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Приводится описание математической модели и результаты численного расчета доли минеральной части угля, улавливаемой в 
топке, применительно к технологии факельного сжигания с удалением шлака в расплавленном состоянии. Рассмотрено влияние 
некоторых конструктивных и режимных факторов на эффективность шлакоулавливания. 


Введение 

Прогнозирование поведения минеральной час- 
ти углей в трактах котельных агрегатов является од- 
ной из актуальных задач для теплоэнергетики, так 
как с этим связано обеспечение надежной работы 
котлов по условиям предотвращения эолового из- 
носа или отложений на тепловоспринимающих по- 
верхностях, поддержания маневренности и несе- 
ния нагрузки, оптимальных условий для процессов 
золоочистки дымовых газов. 

При изучении процессов в котле, определяемых 
свойствами минеральной части угля, моделирова- 


ние условий шлакообразования и компьютерный 
расчет коэффициента шлакоулавливания являются 
ключевыми, поскольку от результата исследования 
зависит надежность моделирования взаимодей- 
ствия минеральной составляющей топлива с пове- 
рхностями нагрева в топке и газовом тракте. Учи- 
тывая определяющую роль аэродинамики для орга- 
низации пылеугольного факельного сжигания в 
энергетических котлах, необходимо исследовать 
шлакообразование в комплексе с моделированием 
двухфазных течений, тепло- и массообмена и горе- 
ния. Для этого применен пакет прикладных прог- 
рамм РІКЕ ЗБ [1], развитый на основе програм- 
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много комплекса СНАІР [2] и использующий в ка- 
честве условия улавливания [3] столкновение час- 
тицы со стенкой, покрытой слоем расплава. 


Математическая модель 

В основу модели положен метод, совмещаю- 
щий эйлеров и лагранжев подходы для описания 
движения и теплообмена реагирующих газов и 
взвешенных частиц в топке парового котла [2] . При 
этом общие уравнения движения, теплообмена и 
горения в газовой фазе описаны на основе эйлеро- 
ва способа представления, т.е. используются прост- 
ранственные уравнения баланса массы, импульса, 
концентраций газовых компонентов и энергии для 
газовой смеси. Лагранжев подход применяется для 
описания движения и тепломассообмена одиноч- 
ных частиц топлива вдоль их траекторий с учетом 
обратного влияния дисперсной фазы на несущую 
среду. Турбулентные характеристики газа рассчи- 
тываются с использованием двухпараметрической 
"к-е модели турбулентности, также учитывающей 
влияние движущихся частиц в несущем потоке. Ра- 
диационный теплообмен представляется в рамках 
Р1 -приближения метода сферических гармоник, 
который показывает хорошие результаты примене- 
ния к пылеугольным топкам [2]. 

Полидисперсность частиц топлива учитывается 
путем выделения основных фракций по функции 
Розина- Раммлера, определяющей рассевочные ха- 
рактеристики угольного топлива после помола. 
Поскольку пылеугольное топливо перед сжигани- 
ем подвергается сушке, считается, что в топку пос- 
тупают частицы, не содержащие влаги. 

Перемещаясь по топочному объему, угольные 
частицы нагреваются за счет радиационно-конвек- 
тивного теплообмена, начинается выход летучих 
компонентов, их воспламенение и горение, а также 
догорание коксового остатка. Предполагается, что 
газовая среда в топке может состоять из химически 
инертных молекулярного азота N 2 , двуокиси угле- 
рода С0 2 , паров воды Н 2 0, а также реагирующих 
0 2 , СО и летучих. 

В этом случае математическая модель включает 
следующие блоки [4]. 


1) Газовая фаза 
Уравнение неразрывности 
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Уравнение энергии 
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Уравнение состояния 
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Здесь р - плотность газовой смеси; Щі= 1,2,3) - 
компоненты скорости газа; х, (/=1,2,3) - декартовы 
координаты; ] щг , /,*,„ - массовые скорости выхода 
летучих при пиролизе угольного топлива и догора- 
ния образовавшегося коксового остатка (гетеро- 
генные реакции); по повторяющимся индексам 
производится суммирование от 1 до 3 (соглашение 
Эйнштейна); С і - концентрации газовых компо- 
нентов; р„ 8с, - турбулентная вязкость газа и турбу- 
лентное число Шмидта (8с, =0,9); ^ ѵ „^ - массовая 
скорость горения летучих, зависящая от скорости 
химической реакции горения углеводородов 
С м Нм-/ ѵЧс * и интенсивности турбулентного смеше- 
ния реагирующих компонентов Д,„, - количе- 
ство 0 2 , необходимое для сгорания 1 кг летучих; к, 
в - энергия турбулентности газа и скорость ее дис- 
сипации; /со - массовая скорость реакции горения 
С0+1/20 2 =С0 2 ; - компоненты вектора радиа- 

ционного теплового потока; Ѳ - интенсивность 
теплообмена между несущей средой и частицами. 
В уравнении энергии моделируется перенос тепла 
за счет конвекции, термического излучения, пос- 
тупления тепла вследствие горения летучих и СО, а 
также обмен теплом с дисперсной фазой; М ] - мо- 
лекулярные веса компонентов газовой смеси; Я, - 
универсальная газовая постоянная. 


2) Дисперсная фаза 

Моделирование движения частиц ведется с ис- 
пользованием лагранжева подхода. При этом весь 
спектр пылеугольных частиц, поступающих в топоч- 
ную камеру, разбивается на А групп по размерам на 
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основе зависимости Розина- Раммлера (или вручную 
на основе данных ситового анализа фракций уголь- 
ной пыли после мельниц). При этом считается, что 
поведение каждой группы частиц может быть оха- 
рактеризовано поведением ее представителя - проб- 
ной одиночной частицы. Таким образом свойства, 
которыми в настоящий момент времени обладает 
частица-маркер, распространяются на все частицы 
группы, которую представляет частица-маркер. Для 
расчета параметров, характеризующих состояние 
частицы-маркера, используется следующая система 
обыкновенных дифференциальных уравнений. 

Уравнения движения частицы-маркера /-ой 
группы [4] 
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где II у - компоненты скорости частиц; Ц=Ц+Щ - 
компоненты актуальной скорости газа; р р - плот- 
ность угля; сі р і — размер частиц /-ой группы; с в . - ко- 
эффициент сопротивления частиц; г, - время их 
релаксации; турбулентные пульсации скорости га- 
за И ) выбираются случайным образом на основе 
принятого нормального распределения Гаусса со 
среднеквадратичным отклонением, равным 2/3 к, 

т.е. |<ЦЗ|=1 к - энергия турбулентности. 

Значение § изменяется с помощью датчика слу- 
чайных чисел в процессе движения частицы в ис- 
следуемом пространстве через отрезок времени 
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булентности. 

Компоненты скорости частиц, не отмеченные 
волной, являются средними скоростями, Ц - турбу- 
лентными пульсациями. Эти и другие параметры по- 
лета частиц-маркеров осредняются после расчета 
достаточного количества траекторий из одной точки. 

Уравнения тепломассообмена и горения частиц 
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Здесь т,- масса одиночной частицы с диаметром 
сІ рі ; с р - теплоемкость угля; () скар - теплота сгорания 
кокса; 0 шр - тепло, необходимое для выхода летучих; 
е р - степень черноты частицы; а сот - коэффициент 
конвективной теплоотдачи одиночной частицы 
(а от 4 і )/Я=2+0,459Ке?- 55 РГ , - я ); Х=р{с/ Рг) - коэффи- 
циент теплопроводности газа; 7) - температура 
/-ой частицы угля. 


При определении массовой скорости выхода 
летучих Мр р используется соотношение (закон Ар- 
рениуса), 


М'7 = к 
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где т ш - начальная масса сухой частицы с диамет- 
ром 4»; / и - массовая доля летучих в сухом угле; 
К а р - предэкспонент; Е шр - энергия активации. 

Скорость горения коксового остатка зависит от 
кинетической скорости химической реакции горе- 
ния углерода и от интенсивности диффузии окис- 
лителя к поверхности частицы, 
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Для К в используется зависимость 
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3) Влияние дисперсной фазы на газовую 

В приведенной выше математической модели 
необходимо определить члены в уравнениях пере- 
носа, которые представляют влияние частиц на ха- 
рактеристики несущей среды: І Ааг , /у, Ѳ, 5 к , 8 Н - 

В использованном здесь смешанном эйлерово-лаг- 
ранжевом способе моделирования аэродинамики, 
тепломассообмена и горения в камере сгорания 
котлоагрегата учет обратного влияния частиц на 
газовую фазу выполняется следующим образом. 

Поскольку решение уравнений будет осущес- 
твляться численно, то рекомендуется использовать 
метод контрольного объема, согласно которому об- 
ласть исследования разбивается на конечное число 
непересекающихся объемов. Значения искомых 
функций внутри каждого объема могут считаться 
однородными. При последовательном расчете тра- 
екторий движения частиц-представителей каждой 
группы с учетом изменения их массы, скорости, 
положения и температуры для каждого контроль- 
ного объема вычисляются значения ] уар , ] Лш , /у, Ѳ, 
8 Ь 5 я по следующим формулам [4]. 
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Здесь 8 Н - поглощение излучения топливными 
частицами, 8 к - дополнительная диссипация турбу- 
лентной энергии, обусловленная влиянием диспе- 
рсной фазы, N -число частиц-маркеров, прошед- 
ших через рассматриваемый контрольный объем, 
имеющий размер К; Ат=Ат!' ар +Ат[' шг - изменение 
массы /-ой частицы в процессе выхода летучих и 
горения коксового остатка за время прохождения 
ею данного контрольного объема; т оі - начальная 
масса частицы-маркера; т ы и т , - начальный и те- 
кущий массовые расходы частиц данной фракции; 
іі - время прохода /-ой частицы через контрольный 
объем V; от, - средняя масса /-ой частицы в конт- 
рольном объеме; А 7] - изменение температуры час- 
тицы в контрольном объеме V. 
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Таким образом, порядок расчета температуры 
шлаковой пленки следующий. Начиная с верхних 
сечений топки, где не образуется шлаковая пленка, 
из (2) получаем <5д= 0 (т.к. 5Ж=0). Если же осаждение 
частиц есть, то расчет производится по вышеприве- 
денным формулам. Заметим, что формулу (3) легко 
можно распространить на случай наклонной пове- 
рхности, для чего следует заменить # на ,§хіпа, где а 
- острый угол между горизонтом и поверхностью. 


4) Шлаковая пленка 


Учет образования и движения шлаковой пленки 
осуществляется с использованием решения Мар- 
шака [5], распространенного на случай простран- 
ственных топочных процессов. Согласно предло- 
женному подходу температура поверхности шлако- 
вой пленки определяется из уравнения теплового 
потока на поверхности 

Чп = ОЕ Т (Т А Ф - Т 4 П ) + а (Т ф - Т п ), 


где а - постоянная Стефана-Больцмана, е т - сте- 
пень черноты топки, Т ф - локальная температура в 
топке, а - коэффициент теплоотдачи, Т п - темпе- 
ратура поверхности движущейся шлаковой пленки. 


Решение задачи об определении средней по се- 
чению пленки скорости имеет следующий вид [5] 
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здесь р ш - плотность шлаковой пленки, - относи- 
тельная вязкость шлака на границе перехода в плас- 
тичное состояние, /и п - вязкость шлака при 

температуре Тр. ц п = РІ т п ) = По ' и ( г > г = - ка- 

сательное напряжение на границе шлаковой пленки. 

С другой стороны, уравнение изменения массы 
шлаковой пленки на участке " ІІр-Боѵт" (верх-низ 
ячейки) в интегральной форме записывается: 

(Рш^ х ^п)оти ~(Рш * Ѵх ^П \ір — 8\Ѵ, ( 3 ) 

где 8\Ѵ - интенсивность осаждения шлаковых час- 
тиц на стенку или в пленку из топочного простран- 
ства на рассматриваемом участке поверхности топ- 
ки. Из (1) и (2) получим: 


Результаты вычислительных экспериментов 

Расчеты выполнены для топки с жидким шла- 
коудалением котла БКЗ-220-100ЖШ, по которой 
имеется обширный банк данных, полученных в хо- 
де натурных экспериментов и физического моде- 
лирования аэродинамики [6]. Топка открытого ти- 
па с прямым вдуванием пыли имеет четыре плос- 
коструйные горелки, установленные на фронтовой 
стене с наклоном к горизонту и разворотом к про- 
дольной оси котла (на 15° для крайних и на 7° для 
средних горелок). За счет взаимной ориентации го- 
релок и стен в нижней части топки создается свое- 
образное взаимодействие плоских струй, приводя- 
щее к возникновению в зоне активного горения го- 
ризонтальной вихревой структуры факела, назван- 
ной разработчиками "аэрошатровой". 

Исходя из вариативности воспроизведения ус- 
ловий работы топочной камеры, вычислительная 
оценка осаждения шлака выполнена несколькими 
этапами: 1) исследовано влияние перераспределе- 
ния угольной пыли по горелкам, а также измене- 
ние вертикального угла оси горелок; 2) проведен 
анализ роли дисперсного состава топлива; 3) ис- 
следовано влияние изменения геометрии топки пу- 
тем установки дополнительного пережима разной 
глубины напротив верхнего среза горелок. При 
этом базовым для сравнения вариантов принят ре- 
жим несения полной нагрузки при работе всех че- 
тырех горелок (сочетание 1 -2-3-4), наклоне горе- 
лок на 15° и при одноплоскостном исполнении 
противоположной горелкам стены топки. Оценка 
коэффициента шлакоулавливания произведена по 
осаждению частиц в пределах футерованных экра- 
нов, ограничивающих объем камеры сгорания. 

Условия подачи аэросмеси в топку оказывают 
значительное влияние на процессы горения, распре- 
деления тепла по теплообменным поверхностям, 
осаждения недогоревших частиц и их золовых остат- 
ков на ограждающих поверхностях топочной каме- 
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ры. Трехмерная визуализация распределения шлако- 
осаждения по стенам при различных вариантах соче- 
тания работающих горелок показана на рис. 1. 

Можно видеть, что наибольшая локализация 
осаждения шлака в камере сгорания достигается при 
рассредоточенном вводе угольной пыли в топку (в ра- 


боте 3 или 4 горелки). В базовом режиме (4 горелки) 
расчетный коэффициент шлакоулавливания состав- 
ляет 0,258. При обеспечении нагрузки котла тремя 
горелками при варианте с двумя крайними и одной 
средней (сочетание 1-2-4) коэффициент шлакоулав- 
ливания уменьшается до 0,245 за счет увеличения вы- 
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Рис. 1 . Визуализация распределения шлакоосаждения в топке котла БКЗ-220- 100 ЖШ при разных сочетаниях задействованных 
горелок 
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носа частиц из камеры сгорания через пространство коэффициент шлакоулавливания возрастает более 
вблизи неработающей горелки. Вариант с большей чем в два раза при работе с одиночной горелкой или 
асимметрией (сочетание 1-2-3) увеличивает коэффи- с двумя центральными горелками (сочетание 2-3). 
циент шлакоулавливания до 0,346 за счет концентри- Влияние угла наклона горелок вниз от горизон- 

рования области аэродинамического взаимодействия тали в визуализированном представлении резуль- 
факела с тыльным экраном с соответствующим уве- татов сказывается на сокращении зоны осаждения 
личением толщины шлаковой пленки. и смещении ее в подовую часть топки (рис. 2). 

Варианты визуализации с двумя и одной горел- 3 целом это сопровождается тенденцией увели- 
кой (сочетания 1, 2, 1-2, 1-3, 1-4, 2-3) являются ги- че ния коэффициента шлакоулавливания (рис. 3, а). 
потетическими с точки зрения реальности несения Однако в области угла наклона 20° и более происхо- 
котлом номинальной нагрузки, однако информа- дих повышение сепарации недогоревших коксовых 
тивны для сравнительного анализа режимов рабо- частиц на под топки с внедрением их в шлаковую 
ты с неполным составом углеразмольных мельниц пленку. На практике это приводит к развитию вос- 
и для выявления участков открытых экранов, по- становления железа в расплаве шлака и "металлиза- 
тенциально опасных по условию шлакования. ции" П ода с угрозой взрывоопасности системы ме- 

Наименыний коэффициент шлакоулавливания ханизированного удаления шлака. Уменьшение уг- 
имеет режим работы с крайними горелками (сочета- ла наклона до 10° также несколько увеличивает ко- 
ние 1-4), однако за счет удлинения траектории факе- эффициент шлакоулавливания за счет появления 

ла обеспечивается более полное выгорание топлива. эффекта "наброса" на противоположный экран при 
По сравнению с этим при прочих равных условиях той же дальнобойности горелочных струй. 
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Рис. 2. Визуализация распределения шлакоосаждения в топке котла БКЗ-220- 100 ЖШ при разных углах наклона горелок к го- 
ризонту 
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Угол наклона горелок, град. 



Тонкость помола В 90 , % 



Рис. 3 . Изменение коэффициента шлакоулавливания ( %) от 
конструктивных и режимных факторов топочного про- 
цесса: а) от угла наклона горелок; 6) от тонкости помо- 
ла угольной пыли (по остатку на сите Р90); в) от глуби- 
ны пережима на противоположном горелкам экране 

Диапазон управляющего воздействия на шла- 
коулавливание за счет изменения тонкости помола 
угля в эксплуатационном диапазоне работы угле- 
размольных мельниц лежит между значениями ко- 


эффициента шлакоулавливания от 0,235 до 0,304 
(рис. 3, б). 

В связи с обсуждением целесообразности суще- 
ственного изменения принципов организации аэ- 
родинамики топки путем реконструкции противо- 
положного горелкам экрана с обустройством и на- 
хождением оптимальной глубины аэродинамичес- 
кого пережима [7] показано (рис. 3, в), что увеличе- 
ние его размера приводит к уменьшению коэффи- 
циента шлакоулавливания вдвое (от 0,258 до 0,125). 

Лишь при перекрытии сечения топки пережи- 
мом на 30 % намечается последующее повышение 
коэффициента шлакоулавливания до 0,147. Выяв- 
ленные эффекты связаны со сложной аэродинами- 
кой, формирующейся вблизи пережима при раз- 
личных его размерах. В частности, дополнительный 
пережим разбивает горелочные струи на два потока, 
один из которых уходит вверх над пережимом, дру- 
гой - в нижнюю часть топки. При увеличении глу- 
бины пережима над ним вплоть до выходного окна 
топки отчетливо проявляется обширная рецирку- 
ляционная зона, осаждение частиц на стенки каме- 
ры сгорания уменьшается, зона повышенных тем- 
ператур (зона горения) смещается к низу топки, 
наблюдается рост скорости газодисперсного потока 
над подом топки. Установка дополнительного пере- 
жима уменьшает общий объем зоны реагирования 
горючей смеси, но увеличивает теплообменную по- 
верхность, и при меняющейся аэродинамической 
картине это приводит лишь к незначительному 
уменьшению температуры на выходе из топки. 

Заключение 

На основании выполненного численного ана- 
лиза коэффициента шлакоулавливания примени- 
тельно к открытой топке с жидким шлакоудалени- 
ем, прямым вдуванием угольной пыли и фронталь- 
ным расположением плоскофакельных горелок 
показаны диапазоны управляющего воздействия 
на осаждение шлака в камере сгорания таких ре- 
жимных факторов как сочетания включенных в ра- 
боту горелок и тонкость помола угольной пыли, а 
также некоторых конструктивных характеристик 
камеры сгорания - угла наклона горелок и глубины 
пережима противоположного горелкам экрана. 
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